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Abstract 

The Brill-Hartle gravitational geon construct as a spherical shell of small am- 
plitude, high frequency gravitational waves is reviewed and critically analyzed. 
The Regge-Wheeler formalism is used to represent gravitational wave perturba- 
tions of the spherical background as a superposition of tensor spherical harmonics 
and an attempt is made to build a non-singular solution to meet the requirements 
of a gravitational geon. High-frequency waves are seen to be a necessary condi- 
tion for the geon and the field equations are decomposed accordingly. It is shown 
that this leads to the impossibility of forming a spherical gravitational geon. The 
attempted constructs of gravitational and electromagnetic geons are contrasted. 
The spherical shell in the proposed Brill-Hartle geon does not meet the regularity 
conditions required for a non-singular source and hence cannot be regarded as an 
adequate geon construct. Since it is the high frequency attribute which is the es- 
sential cause of the geon non-viability, it is argued that a geon with less symmetry 
is an unlikely prospect. The broader implications of the result are discussed with 
particular reference to the problem of gravitational energy. 



To appear in Int. J. Mod. Phys. D 



1 Introduction 



Forty years ago, the geon concept was introduced [jT|: zero rest mass field concentrations 
held together for long periods of time by their gravitational attraction. Such constructs 
were motivated by studies of the motion of bodies in general relativity. More recent 
interest arises from the study of the entropy of radiation and from the analogy between 
electromagnetic geons and quark stars p|. Electromagnetic, neutrino and mixed type 
geons were studied and it was suggested that it should be possible to construct 

a geon from gravitational waves |TD|. Brill and Hartle |TT| (henceforth referred to as 



BH) attempted the construction of a gravitational geon model in detail. Later papers 
(P^^ |T3| - see also assumed the correctness of the BH model. In their approach, BH 
considered a strongly curved static or quasi-static "background geometry" 7^,^ on top of 
which a small ripple h/j^^, resided, satisfying a linear wave equation. The wave frequency 
was assumed to be so high as to create a sufficiently large effective energy density which 
served as the source of the background 7^,^, taken to be spherically symmetric on a time 
average. For their analysis, they took the Regge-Wheeler |]10| (henceforth referred to as 
RW) decomposition of hfj^iy in a spherical background in terms of waves characterized by 
the usual quantum numbers /, m related to the angular momentum operators, and by the 
frequency u. They claimed to have found a solution with a flat-space spherical interior, a 
Schwarzschild exterior and a thin shell separation meant to be created by high frequency 
gravitational waves. With the mass M identified from the exterior metric, there would 
follow an unambiguous realization of the gravitational geon as described above. 

To be complete, however, two conditions must be satisfied. Firstly, the gravita- 
tional geon must be a non-singular solution of the Einstein equations in vacuum. Any 
singularities present would indicate the presence of non-gravitational sources T^^, com- 
pactified into points, curves or surfaces, negating the desired non-singular purely field 
structure. Secondly, the consistency of the solution must be demonstrated, namely that 
the background 7^1, is consistent with the time-averaged effective density constructed 
from hf^i^ as source in the region of non-vanishing /i^,^. Regarding the first condition, it 
is straightforward to show that the junction conditions for regularity are not satisfied 
by the BH solution and hence as it stands, cannot be taken as singularity-free. With 
the first condition violated, there is no basis for proceeding with a consideration of the 
second. 

One might reasonably argue that while the given structure is inadequate as it stands, 
an expansion of the shell region into one of finite extent would reveal a well-posed geon 
solution with both regularity and consistency. Our analysis is sufficiently general to 
include this geometry in which the gravitational field decays sufficiently rapidly at spatial 
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infinity, and to consider also the possibility of geons "leaking" radiation to the exterior. 
Odd high frequency modes in the RW formalism were analyzed in conjunction with a 
static and a time-dependent spherically symmetric background metric 7^,^. It was found 
that the Einstein equations do not allow a solution with the required characteristics 
and hence a spherical gravitational geon cannot exist. While the even mode case or a 
case with a more general geometry than spherical was not yet analyzed, it would be 
unexpected that such a geon could be found when the most primitive case is excluded. 
Moreover, the key factor which leads to the non-existence of the spherical geon is not 
the spatial symmetry but rather the high frequency. This fortifies the expectation that 
the result is general. 



A brief preliminary description of this work was published in [|l^ . The present paper 



provides details of the calculations and an expanded study of the gravitational geon 
problem. It also goes beyond the insufficient generality that was analyzed in the earlier 
paperQ In Sec. 2, we review the basic mathematical formalism for the construction 
of gravitational geons. This is used in Sec. 3 to analyze the proposed BH solution. 
In Sec. 4, we attempt the construction of a non-singular solution for a gravitational 
geon with spherical symmetry and contrast the results with the electromagnetic case. 
It is demonstrated that the Einstein equations do not permit the realization of the 
gravitational geon. In Sec. 5, the results of the previous section are discussed and the 
viability of the proposed electromagnetic geon is critically examined. We conclude with 
a discussion of the potential ramifications of these results with reference to the problem 
of gravitational energy. 



2 Gravitational geons 

We consider the spacetime metric given by 

g/iu = liMu + hf,„ , (2.1) 

where we assume that Qf^i, is asymptotically fiat, that 'y^u is a static, spherically symmet- 
ric, asymptotically fiat metric and hf^i, are small perturbations {\hf^u\ << 1) representing 

"'^We are grateful to to Dr. Paul Anderson for very helpful discussions in this regard. 

^The metric signature is — h ++. We use units in which G ~ c — 1. Greek indices run from to 3 
and Latin indices run from 1 to 3 (apart from Appendix B, where they assume the values 0, 2 and 3). 
A comma and a semicolon denote, respectively, ordinary and covariant differentiation with respect to 
the background metric. The Ricci tensor is given by i?^^. = T^^ ,, - T^^^^ + T^^r;^^ - T'^^r^^. 
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gravitational waves. In a system of Schwarzschild-like coordinates {x"} = {t,r,6,(p}, 
the background metric is given by 



where 
and 



(7^1.) = diag (-e'', e^, sin^ 9 
A = A(r) , u = z/(r) 



(2.2) 

(2.3) 
(2.4) 



Following BH, we represent the most general gravitational wave perturbation h^^y of the 
spherical background as a superposition of tensor spherical harmonics: 



h 



+00 +1 

1=0 m=-l 



dujh^/^';;^\r,9,ifi)e'^' + c.c. 



(2.5) 



This is justified by the fact that the dynamics of the gravitational waves in the present 
context are governed by the linearized Einstein equations around the background 7^,^ 
and therefore a superposition principle holds. Due to linearity, we can restrict ourselves 
to a study of the evolution of the single tensor spherical modes. For ease of comparison 
with the BH paper, we will use the RW set of tensor spherical harmonics (|jlO|; |1T6|-||18||; 
see |TP| for a review and for relations with other sets of tensor spherical harmonics). 
An "even mode" (also called "polar mode" by other authors [Q) in the RW formalism 
is factorized as the product of functions dependent only on time, radius, and angles 
respectively. The angular part is determined by the numbers / and m related to the 
usual scalar spherical harmonics. The even modes have the form 



/ -e^Hoir) Hi{r) 
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where Y^"^{6, (p) are the usual spherical harmonics 0. These modes have parity (— The 
"odd modes" (in the RW terminology - also called "axial modes" ) are given by 
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J 



(2.7) 

and have parity (— )'+^. We will consider the case of odd modes only because of the 
relative ease in computations. It is unlikely that the even modes would produce a con- 
trary result although it would be useful if a follow-up calculation were to be performed 
to verify this conjecture. 

A gravitational geon is defined as a bounded configuration of gravitational waves 
whose gravity is sufficiently strong to keep them confined on a time scale long compared 
to the characteristic composing wave period. It is required that no matter or fields 
other than the gravitational field be present. Although one may consider the possibility 
of strong gravitational waves, and the definition of gravitational geon allows for this 
possibility, in this paper we will restrict ourselves to the case in which the amplitude 
of gravitational waves is small. This permits us to apply the linearized Einstein theory 
to the propagation of each single wave in the background created by the average action 
of all the waves composing the geon. Furthermore, it is required that the configuration 
represented by the metric 7^j, be stable over a time scale much larger than the typical 
period of its gravitational wave constituents, and that the gravitational field becomes 
asymptotically fiat at spatial infinity. Gravitational geons were introduced on the basis of 
the analogy with electromagnetic and neutrino geons in the RW paper and were studied 
in greater detail by BH. Wheeler's method of building an electromagnetic geon was to 
replace the details of the electromagnetic field by the time average of the components of 
the electromagnetic stress-energy tensor. Upon averaging over many modes of oscillation 
of the electromagnetic field, one obtains a stress-energy tensor, and as a consequence, 
a gravitational field and metric which are spherically symmetric. Any given mode of 
oscillation is taken to propagate in the spherically symmetric gravitational field created 
by the rest of the radiation. The attempt to build a geon resembles the construction. 



■^Strictly speaking, the radial functions in Eqs. (2.6) and (2.7) depend on I and m and should be 
labelled accordingly. However, this would result in a cumbersome notation that is preferably avoided. 
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in other fields of pliysics, of a system with many (almost) identical components, each 
of which introduces a negligible perturbation in the dynamics of the whole system and 
has an evolution governed by the averaged action of all the other components. An 
example of such a system in Newtonian theory is a galaxy described by the potential 
created by the mass distribution of many stars (here we neglect dark matter, and the 
fact that a potential-density pair usually describes only a single component of a galaxy, 
and is adequate only for certain types of galaxies |^2|]). Each star gives a very small 



contribution to this potential and its orbit is determined by the global galactic potential. 
Consistent with this idea, it is required that 

= {9,u) . (2.8) 

We also have 

<'v)-(|^)-(fe)-°. 

where ( ) denotes an average over a time that is much longer than the typical gravita- 
tional wave wavelength A ( "Brill-Hartle average"). A mathematically rigourous treat- 
ment of this concept is contained in the paper by MacCallum and Taub |2^. This idea 
has proved very valuable and the averaging process has been used by many authors after 
BH, and is well defined only if it is assumed that the typical wavelength A is much 
smaller than the space and time scale of variation L of the background metric 7^^, (high 
frequency approximation) pT| : 

e=y<<l. (2.10) 

This assumption provides us with a smallness parameter e to be used as an expansion 
parameter. Following [^, we measure times and lengths in units of L so that A = e. 



We have also 

V = 0(e) , (2.11) 
.^^^0(i). (.12) 



"^The term "typical gravitational wavelength" A may be source of confusion to some readers. Since we 
are decomposing the general wave form into an infinite se t of R egge- Wheeler modes, one may think that 



A represents the wavelength of each mode, and that Eq. ( 2.10 ) is only valid if the geon was composed of 
one and only one mode. However, when one is analyzing a general wave form, it is justifiable to assign 
a single parameter describing the scale of variation of the wave form. In the present context, A is the 



scale over which the wave form varies. Equation (2.10) is easily derived from Eq. (2.23) if one keeps in 



mind that h^,^ a ^ e/A etc. (see |£4[) and that A represents the scale of variation of h^^. 
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O 



0(1] 



0(1 



(2.13) 



(2.14) 



In our notation, 0(l)=0(e°). Equation (|2.11|) is derived in [|2T], |2^, It is to be noted 
that, in the most general case of high frequency gravitational waves on a curved space- 
time, two smallness parameters are involved: the dimensionless amplitude of the waves 
and the ratio X/L. These two parameters coincide in the specific case under considera- 
tion, in which the only source of the background curvature are the gravitational waves. 
One can conceive of situations in which more than one parameter arises from the high 
frequency approximation, and these cases have been considered in the literature (see 
e.g. ||2^). However, in these situations, gravitational waves are not the only source of 
curvature. When gravitational waves are the only source of curvature, as in the gravita- 
tional geon, these multiple parameters reduce to the single parameter e. Equation ( |2.13| ) 
implies that the quantum numbers I and m are of order 0(l/e). 



The Ricci tensor can be expanded in the form [11, 



R.p{g) = (7) + R^^l (7, h) + (7, h) + 
where ( [|ll|, |2l|] and references therein) 



(2.15) 



(2.16) 



and h = h"^ 



R. 



(2) 

•a/3 



1 

2 



;/3 



hpT\a ~\~ {hTp\aP ~\~ ho^p-pj- h-pa;Pp ^T/3;ap) 



+ h/3'^'^ {hra;p — hpa-r) — ( ^''^;p ^ 1 (^ra;/3 + ^T/3;a — ^a/3;7 



(2.17) 



?(2) 



The term _R^°j (7) is the Ricci tensor of the background metric 7^1,, whereas 



R)^p and R)^q are, respectively, the parts of the Ricci tensor linear and quadratic in h 



pu 



and their derivatives. In the absence of high frequency waves (or on a flat background), 
hp_i, and their derivatives are all of order 0(e). In this case the superscripts on the 
expansion terms of Eq. ( p. 15 ) also indicate its order in powers of e. However, in the high 
frequency approximation it is clear that contains terms of order 0(l/e) and 0(1) as 
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well as 0(e) Similarly, _R^^J is comprised of terms of order 0(1), 0(e), etc. Solving 
the vacuum field equations 

R^. (g) = (2.18) 

consistently to any order of approximation requires that we set each order in the expan- 
sion parameter e equal to zero. We express Eqs. ( p.l6| ) and ( |2.17|) as 



(7, h) = \e-'] + 



flU 



+ 



i?S(7,/^) = i?SH + ^SH + --- 



(2.19) 
(2.20) 



where R^J^J [e"] denotes the term of order 0(e") in R^J^J. The first order approximation 



is thus 



. 



(2.21) 



The second order approximation requires that terms of order 0(1) be set equal to zero. 
The field equations to this order are 



R^S (7) + R% 



R^^ 



. 



Performing the Brill-Hartle average on Eq. (|2.22| ), one obtains 



R^S (7) 



(2.22) 



(2.23) 



Note that from Eq. 



and hence 



<^ (7, h) 



. 



(2.24) 



(2.25) 



In Eq. ( |2.23| ) the part of the Ricci tensor quadratic in h^^^ and their derivatives has 
been taken to the right hand side and is seen as an effective source term due to the 
gravitational waves. It is important to note that Eq. ( p.23| ) has the potential to lead to 
the description of a gravitational geon only by virtue of the high frequency approxima- 
tion. Under the assumption that gravitational waves are weak but not of high frequency, 
Eqs. ( |2.12| )-( |2.14| ) would not hold and the two terms in Eq. ( p.23| ) would have different 
orders. = O(e^) could never balance R^lil) = 0(1) in this equation. This would 
prevent a priori the construction of a gravitational geon. This point can be understood 
physically by noting that the effective energy density associated with gravitational waves 
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with amplitude h << 1 and frequency u is roughly proportional to {hu)'^. This quantity 
can be of order unity only ii u 1/h » 1. Therefore, it is clear that the high frequency 
approximation is a necessary condition for geon construction in the present context. 

We shall designate as the "geon problem" , the problem of finding a solution {■y^u, h^y) 
to the Einstein equations ( p.21| ), ( p.22| ) and ( p.23|) with the above mentioned properties 
and satisfying the boundary conditions describing asymptotic flatness 

h^j.y — as r ^ +oo . (2.26) 



3 The BH analysis 

To the authors' knowledge the only explicit attempt at gravitational geon construction 
was that of BH. In this Section we review their pioneering approach to the problem and 
critically analyze their work. 

We follow BH in expressing the gravitational wave perturbations in terms of RW 
tensor spherical harmonics. For the sake of simplicity, as done by BH, we restrict 
ourselves to the case of odd modes with zero angular momentum along the z-axis (i.e. 
m = 0). The last assumption eliminates the (/^-dependence from the h^^ functions and 
considerably simplifies the Einstein equations. This can be seen from Eq. ( p.7| ) and 
from the well-known form of the spherical harmonics that we present in Eqs. (|3.4|) , ( p.5|) 
below. Thus, the metric perturbations are Q 

h^,{t,r,e)=n^,{r)Q\e)e-'^\ (3.1) 

where 

K^.(r) = fto(r) (SISI + + ft,(r) (ijij + SlS]) , (3,2) 

e'(9) = .„«^i;=C».„«^%i^. (3.3) 

Here we use the expression of the spherical harmonics 

r'™(^, ^) = V"^P''"(cos d) (m > 0) , (3.4) 

^For ease of comparison with the BH paper, we use a complex exponential to describe the time- 
dependence of the metric perturbations in Eq. (^^). This notation is adequate as long as linear 
quantities in /i^j, and their derivatives are considered, but clearly it is incorrect when the part of 
the Ricci tensor quadratic in h^u and their derivatives enters the discussion. For future reference, we 
use a function of cos {ujt) and sin [uot) instead of a complex exponential in our calculations of Sec. 4. 



8 



¥^"^{6, if) = (-1)" (f'H)* (m < 0) , (3.5) 



im I 



X] 



where C'"^ are normalization constants. Here * denotes complex conjugation and P' 
are the associated Legendre polynomials (which can be expressed in terms of the Leg- 
endre polynomials P\x)). Using the relation P"^(x) = P\x) we obtain 

r^O(^) = C'°P'(cos^) , (3.6) 

from which Eq. (p.3|) follows ^ 

One can now insert the form ( |3.1| )-( |373D of the metric perturbations into the Einstein 
equations ( 2.18| ), obtaining equations for the unknown functions hQ{r) and hi{r). Simul- 



taneously solving Eqs. ( ^.211 ) and ( |2.23| ) for a pair {'y^u, h^^^) then provides a solution to 
the geon problem. 

The correct order of magnitude of the various terms in the Einstein equations is 
determined by Eqs. (|2.11|) - (|2.14|) . The correct order of magnitude decomposition of 



the Einstein equations is absent in [|rT|. While the high frequency approximation was 
assumed in |ll|, it was not incorporated into the calculations. As a result, the authors 
did not obtain the two different orders 0(l/e) and 0(1) in the Einstein equations, using 
a parameter e arising from the high frequency approximation. This is evident from the 
fact that their final equations (lOa)-(lOc) and (14) contain terms of different orders in 
the high frequency limit. In the remaining part of this Section we will show how the BH 
results can be reproduced and we will comment on their proposed geon model. 

The BH equations can only be reproduced in the absence of high frequency waves. In 
terms of a parameter e related to the weakness of the gravitational waves, Eqs. ( |2.11| )- 
( p.l4| ) must be replaced by 



O (V) = O ( ^ ) = O i-^g^ = 0(e) a,P = 0, 3 . (3.7) 




As a consequence of these equations, the Ricci tensor has the form given by Eq. ( p. 15 
where i?Jfj(7) =0(1), = 0(e) and Pg] = 0{e^). To order 0(1) the Einstein 
equations give the well-known equations for a spherically symmetric, static background 
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Note a misprint in the second of the equations (8) in corr esponding to our Eq. (3^). Also to 



be noted is an inconsistency in the notation therein: the form (3.1)-(3.2) for the metric perturbations 
is assumed in p| , but the number m in the definition of the function 0'™ corresponding to our is 
retained. This is inappropriate since it is clear from Eqs. (8) and (9) in [p] that the intention was to 
set m = 0. Otherwise, the function 0'"' would depend on both 9 and ip, which is not the case, and the 
Einstein equations would be much more complicated. 
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(see e.g. p. 300) with vanishing energy-momentum tensor. As far as the order 

0(e) is concerned, only the (0, 3), (1, 3) and (2, 3) components of the Ricci tensor give 
nontrivial results. These components are 



(1) 

■03 



h 



13 



h" 2^^^ ' 

13 ^ "'03 ~ ^^03 



+ ^ (/i03,22 - /i03,2 cot I 



2h 



13 



'■03 



----- 



^ (/il3,22 - hi3^2 cot I 



d(1) 
^23 



h'n,2 - 2^'i3 cot e + /ii3 (A' - u') cot 6 + 



^13,2 / / 



— e " hm cot 9 — 



03,2 



(3.^ 



(3.9) 



(3.10) 



where a dot and a prime denote differentiation with respect to t and r, respectively. 
We now insert the form of the metric perturbations ( p.l|) -( pl3D into the Einstein equa- 
tions (p.21|) and use the following property of the function 9' (see Appendix A): 



f]2p)l jp,l 



After some manipulations we find [| 
h\ + hi 

V ' 



ILO 



2 _ y _ z/^ 
r ~ 2 2 



+ hi 



l{l + l) 



ioj e "hfi + e 



hi 



2u' 



K + Y - A') 



. 



(3.11) 

, (3.12) 

(3.13) 
(3.14) 



^See Ref. for corrections to the radial equations in Refs. |Tl|. Also note misprints in the BH 
Eq. (11) corresponding to our Eq. ( 3.15| ). One of the coefficients of Q in our Eq. (3.19) differs by a 
factor 1/2 from the corresponding one in BH Eq. (14). The sign of the right hand side of our Eq. (3.18) 
is opposite to that in the corresponding BH equation. 
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Following BH we can now use Eq. ( p.l4| ) to eliminate ho from Eq. ( p.l3| ), obtaining the 
second order differential equation for hi{r): 



h'l+h[ 



A")-/(/ + l)^+o^V-'^ + - 



We introduce the variable Q and the Regge- Wheeler coordinate r* defined by 
In terms of these quantities we have 

1 d{rQ) 

IT'O — — 



iuj dr^, 



(3.15) 

(3.16) 
(3.17) 

(3.18) 



anc 



fQ 



+ 



/(/ + 1) 



Q = 



(3.19) 



This Schrodinger-like equation lends itself to the analogy with the dynamics of waves 
propagating in an effective potential |10|, pJ]j . 

At this point BH proceed with the specification of the background metric 




^ (l-2M/r)"^ 



if r < a 
if r > a 
if r < a 
if r > a 



(3.20) 



(3.21) 



where a = 9M/4 and M is the geon mass. This vacuum solution for the background 
metric implies that the effective energy density due to the gravitational waves vanishes 
for r 7^ a. Since the effective energy is positive semi-definite, Eqs. ( p.2q ), ( p.21|) imply 
that 



h 



for r a . 



(3.22) 



An equation similar to Eq. ( |3.19 ) can be derived for the even modes with to = 
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Conversely, if the condition ( p.22| ) is satisfied, the Birkhoff theorem guarantees that the 
metric is Minkowskian for r < a and the Schwarzschild metric for r > a. 

Therefore, in the BH model, gravitational waves are confined to a spherical shell, 
the thickness of which is exactly zero. Apparently, BH meant to build a geon model in 
which the gravitational waves are trapped in a spherical shell which has a nonvanishing 
thickness which is much smaller than its radius. However, their equations do not allow 
for this possibility. To be complete, we examine the viability of a geon with gravitational 
waves confined to a shell whose thickness is exactly zero. It is easy to see that such a 
model is physically meaningless and that the geon problem becomes mathematically ill- 
defined in this case. In fact, the solutions of the radial equations ( |3.12| )-( ]3.15| ) cannot 
be ordinary functions but must be sought in some space of distributions. In Eq. ( p.l5|) . 



the coefficients proportional to v' — A' and v" — \" are not ordinary functions and have a 
mathematical meaning only if they are regarded as distributions. The first of these two 
quantities can be expressed as 



2M 



-1 




z/' - A' = 4Mr-2 j^i _ __ j Oh {r - a) , (3.23) 
where 

if X < 

(3.24) 

is the Heaviside step function. Clearly, the radial derivative of v' — A' can be taken 
only in a distributional sense. Therefore the solutions of the Einstein equations are dis- 
tributions and their domain is some space of test functions which must be specified in 
such a way that the coefficients and the operations involved in the Einstein equations 
are well defined. There is no indication as to the manner in which this functional space 
should be determined. It seems almost certain that, if a meaningful and unambiguous 
mathematical formulation of the problem can be given, the distributional solutions h^y 
cannot be seen as locally integrable functions, but rather must have properties like a 
Dirac delta with support on r = a. Furthermore, the product of distributions is not 
defined and the Einstein equations involving the part of the Ricci tensor quadratic in 
h^v and its derivatives is mathematically meaningless in this case. This destroys the 
possibility of exploring one of the essential features of a gravitational geon. Moreover, if 
the h^j_y are allowed to be distributions, the whole meaning of the linearization around 
the background 7^1,, the condition \h^i,\ << 1, and the estimates of the different orders 
of magnitude in the Einstein equations, become meaningless. The physical interpreta- 
tion of a distributional metric and Riemann tensor is problematic. To appreciate this. 
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one can consider the much simpler case of a metric which does not satisfy the appro- 
priate junction conditions on a spacehke or timehke hypersurface (this is the case 



of the metric 7^,^ given by Eqs. ( ^.201 ), ( |3.2lD and the timehke hypersurface r = a - see 
Appendix B). As suggested by Israel ||2^, and as can be seen from the computation of 
the Einstein tensor for the spherical metric specified by Eqs. (|3.2(]| ), (|3.21| ), a singular 



hypersurface S (in the sense that the first, or the second fundamental form, or both 
are not continuous at S) is associated with nonvanishing T^,^, a source of stresses. The 
definition of a geon, a structure of pure gravitational waves in the absence of matter, 
excludes the use of a background metric which does not satisfy the proper junction 
conditions. If, in addition, the "perturbations" h^i, are allowed to be distributions, the 
consideration of junction conditions loses its meaning, but the argument shows that 
delta-like sources of stresses are included in the problem. Thus, we exclude the case in 
which gravitational waves are confined to a shell, the thickness of which is exactly zero, 
as physically meaningless, mathematically ill-defined, and nonviable. 

The only possible alternative for a geon model in which gravitational waves are 
confined to a spherical shell is the case in which the shell has a nonzero thickness. 
Apparently, BH meant to consider such a model, although this contradicts some of their 
equations. To be specific, let us consider a shell of radius a and thickness Sa described 
by values of the radial coordinate in the range 

6a 6a , ^ 

a--<r<a + - , (3.25) 

where < 6a << a. In order for the geon to be a distribution of pure gravitational fields 
without matter, we must require that the metric tensor satisfies the appropriate junction 



conditions p8[ at the two timelike hypersurfaces S± = {{t,r,6,(p) : r = a ± 6a/2}. This 
guarantees the absence of a real (as opposed to "effective" , i.e. generated by gravitational 
waves) stress-energy tensor T^,^ representing a matter distribution. In this model, the 
modified BH solution would be 

1/9 if r<a-6a/2 



1 - 2M/r if r>a + 6a/2 



(3.26) 



r 1 if r < a - 6a/2 

= , (3.27) 

[ (1 - 2M/r)^^ if r>a + 6a/2 

6a 6a , , 

hf,y = Q if r<a-— , r>a+ — . (3.28) 
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The form of the background metric 7^,^ inside the spherical shell is not given by BH and 
must be determined by solving simultaneously the Einstein equations to the two lowest 
orders for a pair {'j^u, h^^) f^^]- The proper orders of magnitude did not appear in |T1 



as a consequence of neglecting the high frequency approximation, despite the fact that 
this was introduced at the beginning of the paper in order to define time averages. These 
are the reasons why there is only one set of equations in |Tl| mixing different orders and 
a complete solution to the geon problem is not provided. It is natural to ask if such 
a solution based on a spherical shell of nonvanishing thickness is viable. This question 
will be answered in the next Section. 



4 Resolving the geon problem 

In this Section we study the geon problem assuming the high frequency approximation, 
as required, and we take into account the orders of magnitude accordingly. In what 
follows, we solve the geon problem in the case of a spherically symmetric, static and 
asymptotically flat background 'y^u- We solve simultaneously (using numerical meth- 
ods) the Einstein equations to the two lowest orders for a pair {•y^ujh^i,), taking into 
consideration the boundary conditions (p.26|) . We do not assume that the metric per- 
turbations vanish outside a certain radius, but rather solve the system of equations 
throughout space. We consider odd modes only and specialize to m = 0. This is anal- 
ogous to Wheeler's electromagnetic geon. At the end of this Section, the results will 
be generalized to 7^,^ being a time-dependent, slowly varying, spherically symmetric 
background metric. 

4.1 Odd modes 

To avoid problems which can arise in non-linear equations using a complex exponential 
to describe the time dependence of the metric perturbations, it is advantageous to rewrite 
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the individual modes of the RW spherical harmonics in their real formQ for m = 0: 

( Q /io(r)e'(^)ro(t) \ 

hi{r)Q\9)Ti{t) 

Sym Sym 



^ Sym Sym 







(4.1^ 



where 



d 



Q\e) =C'°sin^— P'(cos( 
da 



(4.2) 

To(t) = cos(u;t + 6) , (4.3) 

Ti(t) = sin(c<jt + 6), 6 = constant . (4.4) 

The phase constant 5 can be set to zero without loss of generality, because the phase 
dependence no longer exists upon time averaging. 

The Ricci tensor is computed using Eq. ( ^.16[ ) which, to the dominant order 0(l/e), 
is simplified to (see Appendix C) 



R 



(1) 

'af3 



(4.5) 

(4.6) 

(4.7) 
(4.8) 

from (a,/5) = (0,3), (1,3) and (2,3) respectively. Eliminating hQ^r) from Eq. ( [4.7|) and 
( |4.8|) yield the second order radial wave equation 



Substituting Eq. (f4.1|) into Eq. ( ^3|) yields the three non-trivial equations 

/i'o'(r) - ujh{r) - Pr-^e^ho{r) = , 

uh'oir) + {fr-^e" - u}^) h{r) = , 
h[{r)+uje^--'hoir) = , 



/i';(r) + eM V - - Ui(r) = 



(4.9) 



^Ref. [|5| used Tq — Ti = cos{ujt), which is not the most general case. The functions Tq and Ti are 
constrained by the (0,3) component of the Einstein equations. 
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for the radial function hi{r). It should be noted that in the high frequency limit, the 
angular momentum quantum number / ^ 1, hence Eq. ( 3.11| ) takes the form 



+ /2 = 



(4.10) 



and is used in obtaining Eqs. ( [4.6|) -( ^l9D (note that cot 6dQ^/d6 in Eq. (|3.11|) is of higher 
order than the retained terms in Eq. ( [4.1(j| ) in the high frequency limit). 

The next step is to determine the order 0(1) equations. Instead of evaluating 
Eq. {^^, we can equivalently evaluate 



By defining the Brill-Hartle space-time averaged stress-energy tensor as 



T-iBH 



1 



Stt 



Stt \ ^'^ 



Eq. ( [4.11|) takes the familiar form 



(4.11) 



(4.12) 



(4.13) 



of an effective stress-energy tensor generating the background metric 7^j,. 

The procedure for finding the average over the angular dependence of T^^ has been 
given by Wheeler (see ref. [|l[, p. 520). The results are directly applicable to the gravi- 
tational geon since the discussion covers general T^^. Hence, we have for the time-angle 
average (denoted ()ta) of summed over all N active modes 



o\ 



(2)0 



sin 9 dO 



2\ 



2/. 



SttT, 



TA 
3\ 



3 / 

/ TA 

SttT, 



(I) / T 

sin 6 dO , 
4)2+^(?)3\sin^ci^ 



(I) 1 / T 



for /i 7^ z/ 



(4.14) 

(4.15) 

(4.16) 
(4.17) 
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Here ()t denotes a time average, and is mode number I of the disturbance under 



discussion. Using the mixed form of Eq. ( 2.17|) in conjunction with Eqs. (4J.) and ( [4 .81) , 
and performing the time and angle averaging (see Appendix D), we obtain 



87rT° 



Z 7 Z —U 

— -tu n^e — e 



-A 



2hih'; - h'{ 



(4.18) 



1 1 



^^2gA 



11 2 1 



h't 



2 ^ 2 



h't + Hi 



(4.19) 



We see from Eq. ( [4.16|) that the (2,2) field equation must be identical to the (3,3) 
equation. Neither are necessary to solve the geon problem since the complete system of 
equations must be self consistent. The latter two equations may be derived from the 
(1,1), (0,0) and wave equations. It is of interest to compare Eqs. ( [4.18| )- p.l9| ) with 
their electromagnetic counterparts. The following equations are given in Ref. |l| and are 
not restricted to the high frequency approximation: 



Nl{l + 1) 



TA 2(2/ + 1) 



\rar j 



R 



(4.20) 



87rT; 



o\ 



0/ 



Ar/(/ + l) 



TA 2(2/ + 1) 



— e 



-A 



jdr j \ r J 



1)(4 



(4.21) 



Here R = R{r) is the electromagnetic counterpart of hi{r). It is evident that the similar 
terms (same differential order) in Tq for the gravitational case have a sign difference. 
For the Tq component, one of the three similar terms has a sign difference. These 
sign differences play an important role in the subsequent solving of the gravitational 
geon problem. It will be shown that the sign differences lead to a negative mass for 
the gravitational geon as the only non-trivial solution to the field equations. It is 
understandable that these sign differences should arise: in the electromagnetic case, the 
source is designed at will (subject to the Maxwell equations) and appears on the right 
hand side of the Einstein equations. By contrast, the gravitational source is artificially 
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constructed from terms on the left hand side of the Einstein equations which are shifted 
to the right hand side. 

The left hand side of Eq. { fTWj ) is 



(4.22) 
(4.23) 

The equations for the proposed gravitational geon have the same scale invariance as those 
for the electromagnetic counterpart. Therefore we can introduce the same dimensionless 
measure of radial coordinate, p = uor. By analogy with the electromagnetic case, we can 
define a new measure of potential 



/(p) = v^cu/ii(r), where = -IN{C^'^ 



(4.24) 



and two metric functions L{p) and Q{p) such that 

e-A = 1 _ 2p-^L{p) , 

e^+^^Q\p), 
1 - 2p-'L{p)\ Q\p) . 
Substituting the above into Eqs. ([4.9|), ( [4.18| )-( p:.23| ) yelds the wave equation 



dp*"^ 



. P , 



2L 



f = 



where 



dp* = (^1 - y ^ dp 



and the two background equations 



dL 
dp* 



(1 - 2L/p) 
2Q 



2^ ■' dp*^ [dp* 



1 fd^fY df d^^f 1 2,2^2/. or / ^ 



( 



df 



(4.25) 

(4.26) 
(4.27) 

(4.28) 
(4.29) 



(4.30) 
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dQ 

dp* 

+ 



:i - 2L/p) 

P-2L 



f + f 



d^fl 
dp* dp 



*3 



-p-W(i 




(4.31) 



As in the electromagnetic case, the above equations permit a further change of scale 
without change of form: 



p = bpi, Q = bQi, p* = pI, L = bLi, / = 62/1 



(4.32) 



This allows Qi{0) = 1 and then the scaling parameter can be found by demanding 
Q{oo) = 1. Note that in Eqs. (|4.30| ) and ([4.31|) (and in Wheeler's electromagnetic 
study), an average over a length larger than the characteristic wavelength of the radial 
function f{p) has notheen performed. This would enormously complicate the system of 
equations ( ^4.28| )-( p:.31| ) and was not believed to be necessary in Ref. [|1|. Our approach 
parallels that of Ref. ^ . 

Since the wave equation for the proposed gravitational geon is identical to the one 
for the electromagnetic counterpart, the same arguments hold for applying another scal- 
ing law valid asymptotically for large / (the high frequency of the gravitational waves 
guarantees large /). By making the transformation 



X 



{p* - I) 



(4.33) 



the entire active region of the proposed geon will be described by a range of x of order 
unity. Wheeler [|l| has provided the expansion of the relevant quantities in inverse powers 
of l^^^. They are 



dp* = l^/^dx , 

Pi = 1 + /'/Vo(x) + • • • , 

Li = /Ao(x) + /2/3Ai(x) + /i/%(a;) + --- , 

Qi = l/kix) + r'/^qiix) + r^/^q2ix) + --- 

h = /^/'0(x) + 0i(x) + ri/302(x) + --- . 



(4.34) 



After substituting Eq. ( [4.34| ) into (|4.28| ), expanding in inverse powers of /^^^ and setting 
the lowest three orders to zero, we find the two algebraic equations 



Ar 



(4.35) 
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Ai = qik^ 



(4.36) 



and the differential equation 



+ j{x)k{x)(f){x) = . 



(4.37) 



Here, we have defined 

j{x) = ^ (-ro - 2q2k' + 2A2 + 3qlk^ + Sk^o) (4.38) 
after using Eqs. ( |4.35|) and (|4.36|) . Repeating this procedure for Eq. ( |4.30|) , we obtain 

(4.40) 



dx 



dx 



. 



The expansion of Eq. (|3l|) yields Eqs. (|09|) , (ICTI) and 



dq2 
dx 



+ 



+ 002 + 0^ + \<llk^^ - ^ro02fc-2 = . 



(4.41) 



(4.42) 



By differentiating Eq. iWM) . utilizing Eqs. ( WM . (HI), (H), (CT)-(l4^) and sub- 
stituting 

dvQ 



— = k 

dx 

(derived from Eqs. ( [4.29|) and ([4.34|) ), we obtain the differential equation 

' 7 , \ 2" 



djjx) 
dx 



3- 



1 



2 



dx 



(4.43) 



(4.44) 
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Solving Eqs. ( |4.37| ), ( [4.39| ) and ( [4.44| ) simultaneously for the three functions j{x) 
and k{x) is sufficient for determining the remaining leading terms in Eq. (f4.34|) . It 
should be noted that since Eqs. ( |4.37| ), ( [4.39| ) and ( [4.44D do not explicitly depend on 
X, the system of equations is autonomous. Hence, if j{x) and k{x) are solutions , 
then so are 0(x + a), j{x + a) and k{x + a) where a is a constant. Thus the gravitational 
geon problem is reduced to finding a solution to the system 



dx"^ 



+ j{x)k{x)(t){x) 



dx 2 ^ 







djjx) 
dx 



3- 



k^ 




(4.45) 
(4.46) 
(4.47) 



with the following properties (which are identical to the conditions for the electromag- 
netic geon problem): 

1. For large negative x: The field factor 0(x) — > and k{x) 1. Under these 
conditions dj{x)/dx = 2 or j{x) = 2x. Choosing the integration constant for 
j{x) = 2x to be zero fixes a and consequently defines the origin of x. This removes 
any ambiguity in the start of the integration process. Thus for large negative x, 
(f>{x) satisfies the equation 

^ = 2x0(x) . (4.48) 
The approximate solution as given by Wheeler is 



(x) = I (-2x) 



^1/4 



exp 



-(-2x)3/2] . 



(4.49) 



2. For large positive x: It is required that 
large negative values. 



0, < k{x) < 1 and j(2;) approach 



The only free parameter is the amplitude A of the wave and this must be chosen so 
that the solution fits the boundary conditions. The non-linearity of the problem makes 
it necessary to integrate the system of equations numerically. The integration is started 
at a; = —4. The initial conditions are as follows: 



0(— 4) = 00 = arbitrary 



(4.50) 
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del) 
dx 



1 



,16 
k{-A) = 1 
j(-4) = -J 



+ V8 



(4.51) 

(4.52) 
(4.53) 



Figure 1 shows the result of the numerical integration using a variable step fourth-fifth 
order Runge-Kutta method. 

In order to understand the significance of the curves plotted in Fig. 1, it would be 
beneficial to briefiy review some of the properties of the proposed electromagnetic geon 
solution. The differential equations for the electromagnetic case are 



dx"^ 



djjx) 
dx 



+ j{x)k{x)(f){x) 
dx ' 








(4.54) 
(4.55) 
(4.56) 



The initial conditions for these equations are given by Eqs. ( [4.50|) -( ^.53|) at x = —4. 
The "active region" of the geon is defined to be the range of p where the square 



bracketed combination of terms in Eq. ( 4.28 ) is positive. In this region, the function 
f{p) has oscillatory behaviour. Where the bracketed terms are negative, the behaviour 
of f{p) is exponential growth or decay. The active region can be identified in the x 
coordinate system as the region where the oscillating factor j[x)k{x) is positive. The 
mass of the geon inside radius p is related to the function k{x) in the following way: 



M(p(x)) = -Ao(a;) = 
with b = l/Qi(oo) = k{oo). This implies that 



2b 



1-k' 



< k{x) < 1 or as a; oo 



(4.57) 



(4.58) 



in order to have a positive total mass. 

An eigenvalue solution for Eqs. ([4.54|) -( p:.56|) is sought for which (j){x) — >• as x — > cxd. 
If the amplitude factor A (which translates into an initial choice of the field factor 0o) 
is slightly higher than the desired eigenvalue, (f){x) reaches a minimum and then rises 
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exponentially and becomes singular at some finite x. For an amplitude slightly less than 
the eigenvalue, 0(x) goes to — oo at a finite x. 

Figures 2 and 3 show that the first eigenvalue (characterized by (f){x) having one 
maxima and no local minima) appears to lie between those amplitude factors A, that 
correspond to initial values of 0o at some point in the range 9.790419489 x 10"^ < 0o < 
9.790419490 x 10~^. The mass factor k{x) give a positive mass throughout the integrable 
region (before the singularity) and appears to have a k{oo) value of approximately 1/3. 
The function j{x) is positive only for a limited range in the neighbourhood of x = 1, 
thus identifying the active region. Qualitatively, these results are similar to Wheeler's 
computations. The only main difference between Wheeler's calculation and the present 
one is that his first eigenvalue lies at some point in the range 1.03000 x 10~^ < (pQ < 
1.03125 X 10~^ and the active region starts at a; = 4.05 and ends at x = 6.02. 

For the proposed gravitational geon, the behaviour of the functions is quite different 
from that of the electromagnetic case. This is evident in the differential equations 
themselves. Comparing Eqs. ( [4.55| ) and ( |4.46| ), there is a sign difference on the right 
hand side of the equations. This is the manifestation of the sign difference identified 
earlier in the comparison of the electromagnetic and gravitational stress-energy tensors. 
Proceeding with the analysis of the equations. Fig. 1 shows the results of the integration 
of 0(x), k{x) and j{x) for the initial value 0o = 4.45 x 10~^. The factor j{x)k{x) becomes 
positive at approximately x = and singular at approximately x = 1. Since j{x)k{x) 
never becomes negative again, it implies that there is no end to the active region and as 
a consequence, (plx) remains oscillatory for large x. Of an even more disturbing nature, 
the function k{x) appears greater than 1 for all x and approaches +oo. The implication 
of this is that the mass of the proposed gravitational geon is negative. As 0o is decreased 
(increased), the singular behaviour of k{x) and j{x) moves to larger (smaller) values of 
X, but k{x) remains greater than 1 and j{x) > once it becomes positive. Apparently, 
the only physical solution would be for 0o = 0, which yields k{x) = 1 and hence a 
zero mass gravitational geon. Thus we conclude that it is not possible to construct a 
gravitational geon and the only physical solution is the trivial solution. 



4.2 The time— dependent and stationary cases 

The previous results can be generalized to the case of a time-dependent, spherically 
symmetric background metric 7^,^(t,r), under the assumption that its time variation 
occurs on a scale much larger than the period of the gravitational waves. In this case 
the high frequency approximation and Eqs. ( p.ll|) -( ^.14| ) remain valid. Equation ( p.2|) 
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still holds, but Eq. ( |2.3| ) is replaced by 

\ = \{t,r), u = u{t,r). (4.59) 

As a consequence of the fact that the estimate of the orders of magnitude in the Ein- 
stein equations does not change, we find in this case the same equations that were 
presented above for the odd modes, and the same conclusions apply. If instead, the 
background metric •y^u {t, r) is allowed to vary on a time scale comparable to the pe- 
riod of the gravitational waves, the high frequency approximation does not hold and a 
gravitational geon cannot be constructed, as explained in Sec. 2. This remains valid for 
any time-dependent background metric j^uitjCc) when symmetries are absent, due to 
the fact that our considerations based on Eq. (|2.23|) do not rely on the assumption of 
spherical symmetry. Apart from this argument, the realization of a geon with a rapidly 
varying background metric 7^,^ is problematic for another reason: If a spherically sym- 
metric background is allowed to vary harmonically with frequency Q comparable to the 
frequency of the gravitational waves, one expects a parametric resonance ||3l| for the 
modes with uj = nQ/2, with n = 1, 2, ■ ■ ■ . The strength of the resonance is a maximum 
for n = 1 and decreases rapidly as n increases. In the limit of a static background, the 
resonance phenomenon disappears. Accordingly, on the basis of studies of perturbations 



of black holes and relativistic stars [2C], it is expected that in the case of a station- 
ary axisymmetric background metric describing a rapidly rotating geon, the resonance 
phenomenon between the perturbations and the background metric occurs. In the gen- 
eral case of a time-dependent and rapidly varying background metric 7^,^ (t, x) without 
symmetries, it is not known how to decompose metric perturbations on a complete set 
playing the role of the tensor spherical harmonics in the spherical case, or even how to 
define frequencies in the strong curvature region. However, if such concepts can be given 
a meaning, it seems reasonable to expect some kind of resonance phenomena between 
the background metric and its gravitational wave perturbations. All these resonance 
phenomena certainly do not contribute to the realization of a stable configuration, but 
rather are associated with instabilities that tend to disrupt the system. 



5 Other approaches to the geon problem 

In the previous sections, we have analyzed the BH construct and the models of grav- 
itational geons conceived by Wheeler. However, one can study different models of a 
gravitational geon and different, independent, approaches to the geon problem, which 
are given in the present Section. 
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A first argument, which provides additional intuitive physical insight, is the following: 
We recall our analogy of Sec. 2 between gravitational waves composing a geon and stars 
composing a galaxy. The high frequency approximation required in the geon case has 
a parallel in the case of a galaxy; it corresponds to the requirement that the individual 
stars have a very high velocity. It is clear that such stars would escape from the galaxy 
and would not be trapped by its potential well. A galaxy cannot be built exclusively 
from such stars in rapid motion. In other words, the system would not satisfy the virial 
theorem and would not be bounded. The difference with the gravitational geon case 
is that while one is not obliged to require that stars have a very high velocity when 
constructing a galactic model, the high frequency approximation is necessary for a geon 
and this, in turn, prevents its realization. 

An independent argument to understand the impossibility of a gravitational geon 
is the following: it is well known that, in the limit of high frequencies, gravitational 
waves obey the laws of geometric optics [^, Spatially closed lightlike geodesies 
exist only inside black holes, which necessitate the existence of singularities. Thus, they 
are necessarily inconsistent with the definition of a geon. The null circular geodesic at 
r = 3M in the Schwarzschild geometry is unstable. It is therefore hard to reconcile high 
frequency gravitational waves with stable trapped graviton trajectories in the absence 
of matter. 

The most intuitive model of a gravitational geon is that of a ball of high frequency 
gravitational radiation, which behaves like a perfect fluid with a radiation equation of 
state. It appears that Wheeler was aware of this possibility, but discarded it as non- 
viable and therefore proceeded to study the more complicated models of Ref. [IH, in which 
the waves do not propagate radially and the radiation is not isotropic. "... one naturally 
recalled that some stars derive their energy almost exclusively from particles; others, 
from a mixture of particles and radiation. The extreme limit of a system deriving its 
mass-energy from radiation alone therefore suggested itself. However, with no matter 
to provide opacity and to dam up the radiation against escape, stability could only 
be maintained by excluding all photon orbits in which the motion is purely radial or 
largely radial ..." However plausible this argument may appear, it relies on the 

theory of the stability of Newtonian stars, and cannot be used for a relativistic fluid ball, 
since in general relativity the fluid pressure contributes to the energy density which may 
bind the fluid and this contribution cannot be neglected in a relativistic star dominated 
by radiation pressure. Therefore, it is important to re-examine the possibility of a 
fluid ball made only of gravitational waves. There have been some papers (see 0, 
^ and references therein) analyzing the properties of self-gravitating electromagnetic 
radiation confined to a spherical box. The radiation is taken to be a perfect fluid with 
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equation of state p = p/3. This model is as applicable to high frequency gravitational 
waves as it is to electromagnetic waves, hence the results of these papers will hold 
true for gravitational radiation. In order to build a geon (both electromagnetic and 
gravitational), the constraint of the spherical box (with reflecting walls) would have to 
be removed. This requirement leads to the impossibility of constructing any type of 
geon using a relativistic perfect fluid model. Weinberg |^ has shown that a highly 
relativistic fluid with p = p/3 can never achieve hydrostatic equilibrium in a finite ball 
through gravity alone. The relativistic equation for hydrostatic equilibrium is 

ip + p)^^r.\i-9ooY^'] ■ (5.1) 



For p = p/3, 

P oc (-^oo)"^^+^^/^^- Since p must vanish outside the fluid, goo would have 
to become singular at its surface. If the surface of the ball is allowed to extend to spatial 
infinity, then from Sokolov's work [^, it is not difficult to establish that the mass of the 



radiation ball diverges (see also Ref. |2J], p. 615, ex. 23.10). 

Another work which casts doubt on the existence of electromagnetic geons is that 
of Gibbons and Stewart |3^. They conclude that the Einstein equations do not permit 
asymptotically flat solutions which are both periodic and empty near infinity. This 
precludes the existence of gravitational geons for at least the periodic case. They also 
suggest that the result may be extended to include the case when there is matter near 
infinity, e.g. electromagnetic or scalar radiation. 



6 Discussion and conclusion 

The results of the previous sections were derived by making use of some particular 
gauge conditions that RW imposed in order to set the metric perturbations in the form 



of Eqs. ( p.6|) and (|2.7| ). However, it is clear from their very nature that our results 
are covariant and gauge-independent, since the solution {•y^^, h^^) = {rj^^, 0) that we 
found has an invariant meaning (for example, the vanishing of the curvature tensor is a 
covariant concept). 

Since a spherically symmetric gravitational geon cannot exist due to the fact that the 
high frequency approximation does not allow a solution with the required characteristics, 
one might ask if it is possible to realize a gravitational geon in a configuration with less 
symmetry. We do not expect that such a geon can be constructed when the most 
primitive case is excluded. The main reason for this belief is that the key factor which 
leads to the non-existence of the spherical geon is not the spatial symmetry but rather 
the high frequency. 
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From a mathematical point of view, the main difference between our approach to 
the geon problem, as compared to that of BH, consists in our explicit use of the high 
frequency approximation in conjunction with solving explicitly for the wave and back- 
ground metric functions in a self-consistent manner. We have already seen in Sec. 2 
that this is necessary for the geon problem to be meaningful. In Sec. 4 it was shown 
that the same approximation prevents the realization of a spherically symmetric geon. 

In his papers on geons. Wheeler [0, ^, 0, || describes electromagnetic and neutrino 
geons as systems which are stable on a long time scale, but not absolutely stable, in the 
sense that they "leak" radiation to the exterior. The rate of the leaking is negligible, so 
that a geon is stable for a long period of time. However a secular instability is introduced, 
which seems unavoidable 0. The BH model of a spherical shell with h^y exactly equal 
to zero outside a certain radius excludes such a possibility, and it could be conjectured 
that this might be the reason why their model is not viable, leaving a possibility open for 
the realization of physically more realistic "leaking" geon^. However, this possibility 
is excluded by our calculations. In fact our boundary conditions ( p.26| ) allow for this 
possibility, which in turn is excluded by our results as well. 

It should be noted that while there is a marked contrast between the characteristics 
of proposed gravitational and electromagnetic geons as evidenced in Figs. 1, 2 and 3, it 
is still unclear that even an electromagnetic geon is a viable entity. This is because the 
electromagnetic geon plots do not display true soliton-like confinement. Rather they 
show a region which, by very high refinement of the eigenvalue, is very near vacuum but 
which is followed by an infinite amount of energy beyond this region. It was contended 
that by continued refinement of the eigenvalue, this infinite energy region would itself be 
pushed out to infinity. Firstly, it is not clear that this really is the case and it is possible 
that infinite refinement of the eigenvalue could lead to a convergence of the infinite 
energy regime at some finite distance from the vacuum region. Such would clearly be 
unacceptable as a model of the desired confined energy concentration. Secondly, even 
if it is the case as was conjectured, that continued refinement of the eigenvalue would 
push the infinite energy regime out to infinity, this would not appear to be the idealized 
soliton-like structure that is being sought. 

Traditionally, the geon was conceived as a structure of small-amplitude high-frequency 
gravitational waves compactified to the point where one could describe the resulting met- 



""^^There is inconsistency in ijll] at this point: in that paper it is required that h^^^ (and therefore Q) 
vanishes outside the spherical sheU. However, the Schrodinger-hke equation that is derived there for Q 



(our Eq. ( 3.19| )) imphes a "leaking" geon, as is stated in |T^. In fact, the function Q has a nonvanishing 



tail for large values of the radius, due to the fact that the effective potential barrier is finite. This effect 
is analogous to the well-known tunnel effect in quantum mechanics. 
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ric as the averaged "background" metric induced by the totahty of the waves plus a small 
perturbation due to the local wave presence. This is what was analyzed in the present 
work. It is natural to consider also waves of "large" amplitude in which case lineariza- 
tion is no longer possible nor is it meaningful to envisage a splitting of the metric as 
before. In fact, to assign a measure to amplitude presupposes a standard for comparison 
and in the present work, the background metric served this role. To speak now of large 
amplitude is to consider waves for which there is no longer a discernible "background" 
and hence no standard for comparison of amplitude measure. This leads to the realm 
of exact solutions. One might ask whether an exact wave-like solution of the Einstein 
equations, singularity-free with localized curvature and asymptotically flat, could exist. 
Existing exact wave-like solutions such as the plane waves of Bondi, Pirani and Robinson 
or the cylindrical waves of Einstein and Rosen are not localized and in the second 
case, are also not singularity-free. While it would appear doubtful that solutions with 
the geon-like properties can exist, to our knowledge they are not ruled out. 

Implicit in the gravitational geon concept is the assumption that the gravitational 
field has some particular essential features shared by other fields. Other fields, even in 
their pure states, carry energy. Energy has a mass equivalent and all masses gravitate. 
Thus, given a sufficient concentration of field energy, one could imagine a gravitated 
concentration into a spherical region with the effective mass displayed unambiguously by 
the coefficient of the 1/r part of the asymptotic static vacuum metric. The gravitational 
geon concept is built upon the assumption that the gravitational field itself, even in its 
pure state, will gravitate and thus have the potential to behave as other concentrations 
of matter or fields. Through the years, various authors such as Isaacson ||2l| have dwelt 
upon the similarities between the gravitational and other fields. For example, Isaacson 
has attempted to establish that there is a basis for considering a certain construct of the 
metric as an energy-momentum tensor of the gravitational field which is as substantial 
as a true energy-momentum tensor. However, this requires averaging and under the 
appropriate limits, his construct merges with the energy-momentum pseudotensor, the 
shortcomings of which epitomize the gravitational energy problem. If the gravitational 
field in its pure form really did have the properties which those authors have ascribed 
to it, then it would seem reasonable to expect that a gravitational geon could, at the 
very least in principle, be constructed. However, given the present results, it is worth 
considering alternative ideas. 

Recently, one of the authors |35, 36] introduced a new hypothesis that gravitational 
energy is localized in regions of non-vanishing energy-momentum tensor. The motiva- 
tion derived from the fact that the traditional means by which physicists have identified 
gravitational energy was through the covariant energy-momentum conservation laws. 
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While those laws were extrapolated to produce energy-momentum pseudotensors, im- 
plying densities and fluxes even in vacuum, the fact is that the laws themselves are 
devoid of content in vacuum, producing the empty identity = 0. Given that there 
is a plethora of possible pseudotensors and, as their name implies, they are not really 
tensors, it was suggested |^ that the root of the ambiguity lies in the extrapolation 
of the conservation laws to regions in which they are without actual content. The hy- 
pothesis goes on to propose that the true expression of the gravitational contribution to 
energy is confined to regions of non- vanishing T^j,. In a sense this is the opposite of the 
Isaacson approach in that rather than being satisfied with a construct which reduces to 
the pseudotensor, the new hypothesis suggests that proper localization is realized when 
the pseudotensor is removed. 

Clearly, the realization of a gravitational geon would negate the new hypothesis as it 
would provide an example of a space totally free of true energy-momentum tensor T^j, 
yet exhibit an unambiguous energy content via its asymptotic metric. While one might 
propose exact plane gravitational wave solutions as counter-examples to the hypothesis, 
it is to be noted that these are unbounded fields with questionable relevance to physical 
situations and more directly, these wave solutions can be expressed in Kerr-Schild form 
for which the pseudotensor vanishes in its entirety ||3^. The gravitational geon is a 
direct challenge to the hypothesis and if the geon cannot exist, the hypothesis has 
passed another test. 
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Appendix A: Derivation of Eq. ( |3.11| ) 



We start from the Legendre equation 



d 

dx 



dP\x) 
dx 



+ l{l + l)P\x) = Q 



and note that 



e'(^) =C'°sin^ 



dP\cos( 



de 



(x^ - l) 



dP\x) 
dx 



(A.l) 



(A.2) 
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where x = cos 6. Using 



and the Legendre equation ( |A.1| ), we find the relations 



^ = -/(/ + l)C'°sin^P'(x) , (A.5) 
du 



d^& 



+ 1)C 



/o 



jP'(x) + (x^ - l) 



dO^ 

Using Eqs. and (|OD in Eq. (Q, Eq. (|3ll| ) follows. 



da; 



(A.6) 



Appendix B: Junction conditions for the BH back- 
ground metric 



We consider the Darmois junction conditions |^8|] for the BH background metric on the 
timelike hypersurface S = {{t,r,6,(p) : r = a} separating the regions of the space- 
time manifold U = {(t,r,9,ip) : r < a}, U = {(t,r,9,ip) : r > a}. = {x"} = 
{t, r, 9, (f} and {M*}j=Q 2 3 = f} coordinate systems in U, U and S, respectively 
(note that, in this Appendix, Latin indices assume the values 0, 2, 3 due to the timelike 
character of S). The unit normal to S is directed along the coordinate basis vector dual 
to dr and has components 

n, = 6y/'. (B.l) 

The metric components ■jf^u in U and •y^u in tJ are given by Eqs. (|2.2| ), (|3.20|) and ( p.21|) . 
The first fundamental form of S has components = 'jij. The second fundamental 
form K^iy = rifj^-iy of any hypersurface r =constant has components 

in coordinates Using the Christoffel symbols of a spherically symmetric metric 

(see e.g. |2^), we obtain the only nonvanishing components 

Koo = -'^e''-'/\ (B.3) 
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K22 



r e"^/^ sin^ 9 . 



(B.4) 
(B.5) 



The Darmois conditions p8| require the continuity of the first and second fundamental 
form across S. The first condition is trivially satisfied, while the second is violated. In 
fact, we have 



a 



lim Kqq = 07^ lim Kqq 

r^a~ r—*a+ 

lim K22 = a ^ lim K22 

r^a" r^a+ 3 

lim i^33 = a sin^ ^ lim i^33 

r— >a"f 



16 
27M 



- sm 



(B.6) 
(B.7) 
(B.8) 



where the BH relation a = 9M/4 was used. 



Appendix C: Dominant order in R^^^ 

The second covariant derivatives appearing in Eq. ( |2.16|) are 



h 



-V^ h a VP h VP h 



Symbolically, we express the various quantities in the last equation as follows: 



r = 7^7 = 0(1) 
7^/1 = 0(1), 
{d^)h = 0(e) , 
hdh = 0{e) , 
rdh = 0(1) , 
{dr)h = 0(e) , 
TTh = 0(e) . 



(C.2) 
(C.3) 
(C.4) 
(C.5) 
(C.6) 
(C.7) 
(C.8) 



By using Eqs. ( C.2|) -( |Q!8| ) in ( P.l|) and then, in conjunction with Eq. ( 2.16 ), Eq. ( |4.1|) 
follows. The quantity (/ipr,a/3 + ^a/3,pr — ^™,/3p — hri3,ap) in Eq- ( [1.1| ) contains terms of 



order 0(l/e) as well as terms of order 0(1). We retain only the former ones in the 
linearized Einstein equations to order 0(l/e). 
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Appendix D: Angle average of in the high fre- 
quency limit 

Equations ( [4.14|) -( ^.16|) includes integrating over the angle ip and dividing by the solid 
angle 47r, thus all that is left is evaluating the 6 integrals (see Ref. [|ll). The 6 dependence 
of Ty, comes in three forms 



sin"2 e (Q\e)) , sin-2 e (Q^e)^ and sin-2 e Q\e)Q\e)^22 . 



where 
and 



B\e) = sin ^ (cos ( 
do 



(D.l) 
(D.2) 
(D.3) 



The exact integrals are evaluated below with the last equality being the value used for 
the high frequency approximation: 



sin'^ 9 [B\9)) sin 9de 

j\iYr'^ e (B\e) sinOde 

\m-'^eB\9)B\e)M sin Ode -- 



2Z(Z + 1) ^ 
" 21 + 1 ^ 

2l^{l + lf 
21 + 1 
2l^{l + l) 
2/ + 1 



The normalization constant for Q\9) is found by requiring 

r Q\e)\in9d9d^ = 1 . 
Jo Jo 



Therefore 



10 



2n 



B\e)) sin Ode 



~ 2^ 



4l\l + l] 



1 -1 



(2/- 1)(2/+ l)(2/ + 3) 



Thus the normalization constant is 



10 



(2/-l)(2/ + l)(2/ + 3) 
87r/2(/ + 1)2 



1/2 



(D.4) 

(D.5) 
(D.6) 

(D.7) 
(D.8) 

(D.9) 
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Fig. 1 Results of the numerical integration for the gravitational geon differential equa- 
tions (Eqs. ( [4.45 )-( ^.47|) ). The integration was performed from x = —4 to x = 35 with 



initial condition 0(— 4) = 0o = 1-0 x 10~^. The active region denoted by the region where 
the factor j{x) k{x) is positive extends from approximately x = to oo. Consequently 
cj){x) oscillates out to oo and the mass of the gravitational geon (proportional to (1 — fc^)) 
is negative for all x. Other values of 00 do not qualitatively change the behaviour of the 
functions j{x) or k{x). 

Fig. 2 Results of the numerical integration for the electromagnetic geon differential 
equations (Eqs. (^34D-(^36D). The initial value of 0(-4) was 0o = 9.790419490 x lO'^. 



The integration started at x = —4 and could not proceed beyond approximately x = 11. 
The active region began at approximately x = 0.1 and ended at x = 0.2. Note that 0(x) 
exhibits singular behaviour at approximately x = 11. 

Fig. 3 Results of the numerical integration for an electromagnetic geon with an initial 
value of 00 = 9.790419489 x 10~^. For this case, 0(x) approaches —oo at approximately 
X = 11. 
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